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H ■ Abstract 

We investigate the acceleration of electrons via their interaction with electro- 
static waves, driven by the relativistic Buneman instability, in a system dominated 
by counter-propagating proton beams. We observe the growth of these waves and 
their subsequent saturation via electron trapping for a range of proton beam ve- 
Q^l locities, from 0.15c to 0.9c. We can report a reduced stability of the electrostatic 

wave (ESW) with increasing non-relativistic beam velocities and an improved wave 
stability for increasing relativistic beam velocities, both in accordance with previous 
I findings. At the highest beam speeds, we find the system to be stable again for a 

Q^' period of ~ 160 plasma periods. Furthermore we observe a, to our knowledge, pre- 

lyj . viously unreported secondary electron acceleration mechanism at low beam speeds. 

We believe that it is the result of parametric couplings to produce high phase ve- 
locity ESWs which then trap electrons, accelerating them to higher energies. This 
(~| '. allows electrons in our simulation study to achieve the injection energy required for 

C^l Fermi acceleration, for beam speeds as low as 0.15c in unmagnetised plasma. 
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1 Introduction 



Supernova remnant (SNR) shocks, created by the interaction between an expanding 
supernova blast shell and the ambient medium, are believed to be a significant source 
of cosmic rays (Reynolds |20nij : Lazendic et al. |2nn4j ). with energies up to 10^^ 
eV, the knee of the cosmic ray spectrum (Nagano & Watson j2UUU| : Aharonian et 
al. [20011; Volk et al. fl988j V 

First order Fermi acceleration has been proposed as a mechanism for the accel- 
eration of electrons to ultra-relativistic velocities in SNR shocks (Bell |1978aj ; Bell 
|1978bj : Blandford & Ostriker jl978j ). Electrons gain energy by repeated crossings 
of the shock front, and by their scattering off MHD waves on either side of the shock. 
The orientation of the ambient magnetic field B relative to the shock normal has 
implications for the efficiency of electron acceleration (Galeev jl984p . We consider 
the case where the ambient field is orthogonal to the shock normal, this geometry 
provides an efficient electron acceleration mechanism for high Mach number shocks 
(Treumann &: Terasawa |2nnij ). First order Fermi acceleration at such shocks re- 
quires a seed population of electrons that have Larmor radii comparable to the shock 
thickness. Since the shock thickness is, for perpendicular shocks, of the order of the 
ion Larmor radius, the electrons of this seed population must have mildly relativistic 
initial speeds. A kinetic energy comparable to 100 keV is believed to be sufficient 
(Treumann & Terasawa "2001 ). Such electrons are unlikely to be present, neither 
in the interstellar medium (ISM) nor in the stellar wind of the progenitor star of the 
supernova, but may be created by a pre-acceleration mechanism at the SNR shock 
front. This pre-acceleration, commonly referred to as the injection problem, is not 
well understood. 

As the shock front of a supernova remnant (SNR) expands into the ISM, it 
reflects a substantial fraction of the ISM ions as observed in simulations (Shimada 
& Hoshino |2000j : Schmitz et al. |2002j ) and in-situ at the Earth's bow shock 
(Eastwood et al. |2005j ) . If the shock normal is quasi-perpendicular to the magnetic 
field in a high Mach number shock, as many as 20% of the ions can be reflected 
(Sckopke et al. |1983j : Galeev |1984j : Lembege & Savoini [T992] : Lembege et al. 
j2()()4j ). The reflected ions form a beam that can reach a peak speed comparable 
to twice the shock speed in the ISM frame of reference as, for example, discussed 
by McClements et al. ( jl997| ). The shock-reflected plasma particles are a source of 
free energy, similar to the shock-generated cosmic rays (Zank et al. jl990j ^ which 
are also thought to heat the inflowing plasma. However, the shock reflected ion 
beam is considerably more dense than the cosmic rays and each particle carries 
less energy. The developing plasma thermalisation mechanisms are thus likely to 
be different. Since binary collisions between charged particles in the dilute ISM 
plasma are negligible, these ion beams relax by their interaction with electrostatic 
waves and electromagnetic waves. In what follows, we focus on the interaction of 
electrons with high-frequency electrostatic waves (ESWs) that are driven by two- 
stream instabilities. 

Recent particle-in-cell (PIC) simulation studies (Shimada & Hoshino |2003j : Shi- 
mada k Hoshino j2()()4j : Dieckmann et al. j2()()()j : Dieckmann et al. 2004 a; : Mc- 
Clements et al. |20Ulj ) have examined these mechanisms with a particular focus 
on how and up to what energies the ESWs driven by non-relativistic or mildly 
relativistic ion beams can accelerate the electrons in the foreshock region. 

The maximum energy the electrons can reach by such wave-particle interactions, 
depends on the life-time of the saturated ESW and on the strength and the orien- 
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Fig. 1: As the SNR shock expands it reflects a fraction of the ISM protons. These protons 
move back into the upstream region and form beam 1. The upstream B rotates beam 
1 which returns as beam 2. The simulation box covers a one dimensional region of x in 
front of the shock sufficiently small that we can assume it to be spatially homogeneous 
and one dimensional. 



tation of B. Initially a stable non-linear wave known as BGK mode (Bernstein et 
al. jl957j : Manfredi jl997j : Brunetti et al. j2UUUj 'l develops if the plasma is unmag- 
netised or if the magnetic field is weak (McClements et al. |2n01j : Dieckmann et 
al. j2UU2j : Eliasson, Dieckmann & Shukla {2005) 1. Such modes are associated with 
phase space holes - islands of trapped electrons. These BGK modes are destabilised 
by the sideband instability, a resonance between the electrons that oscillate in the 
ESW potential and secondary ESWs (Kruer et al. |1969j : Tsunoda & Malmberg 
jl989j : Krasovsky |1994j '). This resonance transfers energy from the trapped elec- 
trons to the secondary ESWs. The initial BGK mode collapses, once these secondary 
ESWs grow to an amplitude that is comparable to that of the initial wave. 

Many previous simulations of two-stream instabilities in the context of electron 
injection and of shocks have employed PIC simulation codes, which suffer from 
high noise levels (Dieckmann et al. |2004cj l and from a dynamical range for the 
plasma phase space distribution that is limited by the number of computational 
particles. It is thus possible that certain instabilities can not develop, due to a 
lack of computational particles in the relevant phase space interval, or that the 
phase space structures (e.g. BGK modes) are destabilised by the noise (Schamel & 
Korn |1996j ) and may benefit from a plasma model based on the direct solution of 
the Vlasov equation. Such effects have been demonstrated in the non-relativistic 
limit, for example by Eliasson, Dieckmann &: Shukla f |2005j 'l and Dieckmann et 
al. ( |2004bj ). in which the results computed by PIC codes have been compared to 
those of Vlasov codes. The particle species in a Vlasov code are represented by a 
continuous distribution function, conventionally evolved on a fixed Eulerian grid, 
rather than by simulation macro-particles. The comparison of results from these 
two methodologies has shown significant differences in the life-time of the BGK 
modes both for unmagnetized and magnetized plasma. 

SNR shocks typically expand into the ISM at speeds ranging between a few and 
twenty percent of c (Kulkarni et al. |1998j b The shock-reflected ions can thus 



reach a speed of Vh ~ 0.4c if the reflection is specular (McClements et al. jl997j ) 
or even higher speeds if we take into account shock surfing acceleration (Ucer & 
Shapiro |2()()lj : Shapiro & Ucer |2003j ) . The proton-beam driven ESWs have a 
phase speed similar to vt (Buneman |1958j : Thode & Sudan jl973j l. Since the 
ISM electrons can be accelerated to a maximum speed well in excess of the phase 
speed of the ESW (Rosenzweig jl988j ^. we must consider relativistic modifications 
of the life-time of the BGK modes in SNR foreshock plasma. It has been found 
by Dieckmann et al. f |2UU4aj ) that the BGK mode is stabilised if the phase speed 
of the ESW is relativistic in the ISM frame of reference. The likely reason is that 
the change in the relativistic electron mass introduces a strong dependence of the 
electron bouncing frequency on the electron speed in the ESW frame of reference. 
This decreases the coherency with which the electrons interact with the secondary 
ESWs and thus the efficiency of the sideband instability. The stabilization is clearly 
visible in PIC simulations for ESWs moving with a phase speed of 0.9c (Dieckmann 
et al. 20Q4aj ). For speeds < 0.7c the relativistic modifications of the BGK mode 
stability have been small in the PIC simulations (Dieckmann et al. |2004aj ) . The 
better representation of the phase space density afforded by a Vlasov simulation 
may, however, yield a different result and it requires a further examination. 

We thus focus in this work on the modelling of electrostatic instabilities by means 
of relativistic Vlasov simulations and we assess the impact of these instabilities as a 
potential pre-acceleration mechanism. To this end, we neglect magnetic field effects 
and consider only the one-dimensional electrostatic system, which is equivalent to 
the approach taken by Dieckmann et al. ( ^2004a] ) but at a much larger dynamical 
range for the plasma phase space distribution. This description serves as simple 
model for the stability of ESWs by which we identify similarities and differences 
between the results provided by PIC and Vlasov simulations. Future work will 
expand the simulations to include magnetic fields, which introduces electron surfing 
acceleration (ESA) (Katsouleas & Dawson jl983j : McClements et al. |.200T] ) and 
stochastic particle orbits (Mohanty & Naik |1998j ). and multiple dimensions. 

More specifically, we examine by relativistic electrostatic Vlasov simulations 
(Arber &; Vann |2002j : Sircombe et al. .2005. ) how the BGK mode life-time and its 
collapse in an unmagnetised plasma depend on vi, and thus on the phase speed of the 
ESW. The purpose is twofold. First, we extend the comparison of results provided 
by PIC and by Vlasov codes beyond the nonrelativistic regime in Dieckmann et al. 
f |2004bj l. We perform Vlasov simulations for initial conditions that are identical to 
those in Dieckmann et al. f 20 04aj ) where a PIC simulation code (Eastwood jl991j ) 
has been used. We find a good agreement of the results of the relativistic Vlasov 
code and the PIC code for relativistic f & and with the corresponding result provided 
by the nonrelativistic Vlasov code (Eliasson ^20(7^ : Dieckmann et al. |2004bj ). We 
thus bring forward further evidence for both, an increasing destabilisation of BGK 
modes for increasing nonrelativistic phase speeds of the wave and a stabilisation 
for increasing relativistic phase speeds. Secondly we want to exploit the much 
higher dynamical range of Vlasov simulations to examine the wave spectrum and the 
particle energy spectrum that we obtain by the considered nonlinear interactions. 
We find that the collapse of the BGK modes couples energy to three families of 
waves. Firstly a continuum of electrostatic waves that move with approximately 
the beam speed. These waves are connected to the turbulent electron phase space 
flow. Secondly we find for high beam speeds the growth of quasi-monochromatic 
modes with frequency comparable to the Doppler-shifted bouncing frequency of the 
trapped electrons in the wave potential. These waves would be sideband modes 



(Kruer et al. |19H9j ). Thirdly we find waves that do not have a clear connection 
to any characteristic particle speed. We believe that these modes are produced 
by parametric instabilities. These modes can reach phase speeds well above the 
maximum speed the initial trapped electron population reaches, and they grow to 
amplitudes at which they can trap electrons, i.e. a BGK mode cascade to high 
speeds develops. By this trapping cascade, the electrons can reach momenta well 
in excess of those reported previously (Dieckmann et al. |2nnnj : Dieckmann et 
al. |2UU4aj ) . This result would imply that the ion beams, that are reflected by 
shocks that expand at speeds comparable to SNR shocks, can accelerate electrons 
to energies in excess of 100 keV, at which they can undergo Fermi acceleration to 
higher energies. 

2 The physical model, the linear instability and the 
simulation setup 

2.1 Physical model 

We consider a small interval of the ISM plasma just ahead of the SNR shock (see 
Fig. ^ and we treat the ISM plasma as an electron proton plasma with a spatially 
homogeneous Maxwellian velocity distribution. We place our simulation box close 
to the SNR shock, so that the shock-reflected ISM protons can cross it. The protons, 
that have just been reflected, constitute beam 1, beam 2 represents the protons that 
return to the shock after they have been rotated by the global foreshock magnetic 
field. Our model is in line with that discussed, for example, by McClements et al. 
f |1997p and solved numerically, e.g. for unmagnetized plasma by Dieckmann et al. 
r |2000l I2()()4a[ I2004bp and for magnetized plasma by McClements et al. f [200T] ^ 
and by Shimada & Hoshino f j2004j ). 

We set -6 = 0, and thus, exclude electromagnetic instabilities, e.g. Whistler 
waves (Kuramitsu & Krasnoselskikh 200.5 ^ ) . MHD waves, and electrostatic waves 
in magnetised plasma, e.g. electron cyclotron waves. However, by this choice we 
decouple the development of competing instabilities and we can consider them sep- 
arately. In this work we focus on ESWs and the nonlinear BGK modes, which are 
important phase space structures in the foreshocks of Solar system plasma shocks 
(Treumann & Terasawa |2001j ). 

The system is, with the choice B = 0, suitable for modelling with an electro- 
static and relativistic Vlasov-Poisson solver. The size of the simulation box is small 
compared to the distance across which the beam parameters change. Thus we can 
take periodic boundary conditions for the simulation and spatially homogeneous 
Maxwellian distributions for both beams. Both beams have the same mean speed 
modulus, \vb\, but move into opposite directions which gives a zero net current in 
the simulation box. We take a higher temperature for beam 2 than for beam 1 to 
reflect the scattering of the beam protons as they move through the foreshock. We 
show the velocity distributions in Fig. [21 

The choice of i? = is a critical limitation of the model. The initial conditions 
described above could not fully describe a perpendicular shock and with B = 
one could not account for the presence of a returning proton beam. We assume 
that, while the field is neglected over the simulation box, the field outside the box 
is sufficient to produce the proton beam structure described. Therefore, the results 
presented here are not directly applicable to the foreshock dynamics of high Mach 
number shocks. However, they provide an overview of processes that require a large 
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Fig. 2: The initial velocity distribution showing the ISM electrons and protons, proton 
beam 1 and proton beam 2. The expanding SNR shock reflects a fraction of the ISM 
protons which form beam 1. The upstream B rotates beam 1 to create beam 2. Since the 
beam has been scattered by ESWs on its path through the foreshock, its thermal spread 
has increased. Both beams can grow ESWs with a speed similar to Vb which saturate by 
their interaction with the ISM electrons. By assuming that the plasma parameters do 
not change over the small box size we consider, we can represent this plasma by the 4 
spatially homogeneous species shown. 

dynamical range for the plasma phase space distribution and which may, therefore, 
not have been observed in previous PIC simulations. Our system is applicable 
to parts of the foreshock of perpendicular shocks, where magnetic field fluctuations 
(Jun & Jones jl999j ) cause the magnetic field to vanish or to be beam aligned. It may 
also apply to the field aligned ion beams that are observed in the foreshock region 
of the Earth's bow shock (Eastwood et al. ^200^) and which may be present also at 
SNR shocks. Here, the second proton beam in our initial conditions takes the role 
of the return current in the plasma which is typically provided by all plasma species 
(Lovelace & Sudan |1971j ). Interactions between two BGK modes in unmagnetized 
plasma affect primarily the velocity interval that is confined by the phase speeds of 
the two waves (Escande jl982| '). Our results, which focus on the developing high 
energy tails of the plasma distribution, may thus not depend on the exact setting 
of the initial return current and may be more universally applicable. 

2.2 Linear instability 

1 /2 

We introduce the plasma frequency of each species i as ujp^i = {e^rii/mieo) where 
e, rii, rrii and eo are the magnitude of the elementary charge, the number density 
of species i in the rest frame of the species, the particle mass of species i and the 
dielectric constant. In what follows we define all ujp^i in the box frame of reference. 




Species 1 and 2 are the ISM electrons and protons. The species 3 is beam 1 and 
species 4 is beam 2. The large inertia of species 2 compared to species 1 impHes, 
that its contribution to the linear dispersion relation can be neglected. The linear 
dispersion relation then becomes in the cold plasma limit 



^ + — 2 + —. — 2 -1 = (1) 

^ Y{vb){uj - Vbk) j^{vh){uj + Vbk) 

During the linear growth phase, each proton beam will grow an ESW by its stream- 
ing relative to the electrons. Both waves are well separated in their phase speed 
and the linear dispersion relation can be solved separately for each ESW by neglect- 
ing either the term with iVp^3 or ojp^4 in Eq. ^ (Buneman jl958j : Thode & Sudan 
|1973j ). Its frequency in the box frame is uJu ~ <^p,i, its wave number is ku ~ ^p,iVb 

1 /3 

and its growth rate is ~ {3\/^L0p /j{vb). We take a number density 

ratio of n^/ni = 0.2'y{vb) and set 714 = 723 and n2 = ni — — in the box 
frame of reference, which is representative for a shock with a high Mach number 
(Galeev |1984j ). This density ratio gives cop^i/cop^^ ^ 96'y{vb)~^'^ and a growth rate 
of w 0.033/7(ffe)^''^. The growth rate is reduced by increasing temperatures of the 
plasma. The wave length of the most unstable ESW is Xu = 27r /ku- The sideband 
instability couples energy to modes with k < ku, for non-relativistic phase speeds 
of the wave (Krasovsky ^94j). We thus set the simulation box length to L = 2Xu 
to resolve more than one unstable sideband mode, which leads to a wave collapse 
(Dieckmann et al. fMK)\ ). 

This L is short compared to the typical size of the foreshock which justifies our 
spatially homogeneous initial conditions and periodic boundary conditions. We use 
the amplitudes of the initial ESW and that of the sideband modes as indicators 
for the wave collapse. We define I as the index of the simulation cell, xi = /Ax, 
kj = 2Trj /NxAx and N^, as the number of simulation cells in x. The index m refers 
to the data time step mAt, where At is the time interval between outputs rather 
than the simulation time step. We Fourier Transform the spatio-temporal ESW 
field E{x,t) as 



E[kj, tr. 



(2) 



Nx~^\'^E{xi,tm) exp (ikjXi) 
1=1 

The amplitude of the ESW with k = ku at the simulation time tm is then given by 

E{k2,tm)- 

To analyse the nonlinear and time dependent processes developing after the 
saturation of the ESW, we introduce a Window Fourier Transform. We define the 
window size, in time steps, as Nt and we introduce uJn = nAuj = n{2'K /NtAt) as 
the frequency defined in a Fourier time window with a size NfAt. This frequency 
is limited by the sampling theorem to — vr/At < tj„ < -k j At. The Window Fourier 
transform can be written as 

W^i(%,*m) = N^~^ E{xi,tm) exp (ikjXi) 
1=1 

s+Nt/2 ^ 

Wi{kc,ti)ex.p{iu;nti) 

l=s-Nt/2+l 



W2{uJn,ts,kc) = Nf 



(3) 



where the Fourier transformed time interval is small compared to the total length 
of the time series, which is IdN^At, and where kc is a fixed wave number. 



2.3 Numerical Simulation 



In the absence of a magnetic field the one dimensional relativistic Vlasov-Poisson 
system of electrons and protons is given by the Vlasov equation for the electron 
distribution function fe{x,p,t) = fi 

'4 + ^f-eEf = 0, (4) 
at rrie^ ox op 

the Vlasov equation for the proton distribution function fp{x,p,t) = /2 + /s + /4 

^ + ^f^+ei.f^ = 0, (5) 
ot rup^ ox op 

and Poisson's equation for the electric field 

The Vlasov-Poisson system is solved using a parallelised version of the code detailed 
in Arber & Vann ( |2UU2j ) and Sircombe et al. ( |2UU5j ). This uses a split Eulerian 
scheme in which the distribution functions (/e, fp) are calculated on a fixed Eulerian 
grid. The solver is split into separate spatial and momentum space updates (Cheng 
&; Knorr 1976 ). These updates are one dimensional, constant velocity advections 
carried out using the piecewise parabolic method (Colella & Woodward jl984j ^. 



2.4 Initial Conditions 

Throughout our simulations we adopt a realistic mass ratio, mp/nie = Mr, of 
Mr = 1836 and resolve L by 512 cells in the x— direction, each with length Ax = 
7rffe/128u;p^i. To accurately resolve the filamented phase space distributions of the 
particles that result from the sideband instability (McClements et al. j2UUlp . we 
use a momentum-space grid with between 4096 and 16384 grid points {Np). This 
ensures that the momentum grid spacing for each species, Ap^^i, is small. Specifi- 
cally, Apx^i < 0.1miVth,i where v^f^- = KTi/nii is the thermal speed of species i and 
where k and Tj are the Boltzmann constant and the temperature of species i. 
We set the temperatures of the four species to Ti = 5.4 x 10^ K, T2 = = lOTi 
and T4 = 100r2. We thus obtain vth,i = 10~^c. Each species is described by a 
Maxwellian momentum distribution of the form 

ft (P) = Ci exp ( 1 (7) 

where po is the initial momentum offset, zero for species 1 and 2, ±Vh/{l — v'^/cP)^'^ 
for species 3 and 4 respectively. For each species the constant Ci is calculated to 
ensure that / fidp = rii and that / (/2 + /s + fi) dp = rii, so there is no net charge 
in the system. 

In order to excite the linear instability, we add to the proton distribution (species 
2, 3 and 4) a density perturbation at the most unstable wavenumber of the form 
n' = a cos kuX, where a is small, typically of the order of 1% of the background 
density. The minimum beam speed of Vb = 0.15c we examine here corresponds to 
the slow beam in Dieckmann et al. ( |2004bj ). The maximum beam speed Vb = 0.9c 
equals the maximum beam speed in Dieckmann et al. ( |2004aj ). 




Fig. 3: The logarithmic amphtude of the most unstable mode E{ku,t) is plotted 
against time, in normalised units, for a variety of beam speeds. The dashed lines on 
each plot represent a linear fit over the region of exponential growth. Using this lin- 
ear fit we obtain growth rates for the relativistic Buneman instability, in normalised 
units, oi n ^ 0.0256, 0.0264, 0.0270, 0.0246, 0.0201 and 0.0137 for beam speeds of Vb = 
0.15c, 0.2c, 0.4c, 0.6c, 0.8c and 0.9c, respectively. 



3 Simulation results 

Simulations use a system of normalised units where time is normalised to uj~\, space 
to cojpl and momenta to "yrrieC. It follows that the electric field is normalised to 
LVp^icnii/e. Where appropriate in the figures, we re-normalise x in units of 2Tr/ku, 
the most unstable wavelength for a given initial beam speed. This ensures that the 
spacial units are identical for all beam speeds. 



3.1 Growth and non- linear Saturation of the Bune- 
man instability 

Figure 01 shows the initial growth stage of the most unstable mode, E(k = ku), 
for a range of initial beam velocities. From these we estimate the growth rate of 
the instability (Q), in normalised units, to be 0.0256,0.0264,0.0270,0.0246,0.0201 
and 0.0137 for beam speeds of Vb = 0.15c, 0.2c, 0.4c, 0.6c, 0.8c and 0.9c, respectively. 
These compare favourably with the linear theory. Writing a{vb) as the ratio between 
observed and theoretical growth rates we find; a(0.15c) 0.79, a(0.2c) 0.81, 
a(0.4c) « 0.87, a(0.6c) ^ 0.87, a(0.8c) « 0.86 and Q!(0.9c) « 0.72. As explained in 
Dieckmann et al. ( j2()04aj ). where a similar systematic reduction has been observed 
in PIC simulations (in this case by 15-20%), this might be connected with the 
fast growth rate of the instability itself since it results in a considerable spread in 
frequency for the unstable wave. This makes the treatment of the instability in 
terms of single frequencies (uJu) inaccurate. In truth, the growth rate should be 
lower since the energy of the unstable wave is spread over damped frequencies. 

The ESWs saturate by the trapping of electrons and the formation of BGK 
modes. This is shown in Fig. 0J for the case of Vb = 0.2c, and Fig. [5] for Vb = 
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Fig. 4: Contour plot of log(/e > 10~^) at time t = 150 for a system with an initial beam 
velocity of Vb = 0.2c. The exponential growth of the ESW is in the process of saturating 
via the trapping of electrons, forming a BGK mode at p ~ +0.2.. The thermal spread 
of beam 2 inhibits its growth slightly, allowing the system to become dominated by one 
ESW. Hence, we do not observe electron trapping around p ~ —0.2. Here x is given in 
units of 2TT/ku- 

0.9c. Here we see the appearance of phase-space holes characteristic of particle 
trapping. While the saturation mechanism is the same in both cases, for the high 
velocity beam we observe the development of two counter-propagating BGK modes. 
This is because at lower beam speeds the ESW instability driven by the cooler 
beam (at +Vb) dominates whereas for higher, relativistic, beam speeds forwards and 
backwards propagating ESWs grow in unison. At lower beam speeds the increased 
temperature of beam 2 (species 4) is more significant, since the thermal velocity of 
the beam is a larger proportion of the beam speed than is the case at Vh = 0.9c. 
Thus, the growth rate for the hot beam is sufficiently reduced for the system to 
be dominated by the growth of the cooler proton beam (species 3). We do observe 
the growth of a counter propagating ESW which begins to trap electrons after t = 
150. However, the final momentum distribution is clearly dominated by electrons 
accelerated by the ESW with positive phase velocity. At Vb = 0.9c, the thermal 
spread of both proton beams is negligible in comparison to the beam velocities and 
we observe the growth and saturation of ESWs associated with both beams. 

3.2 Wave collapse and secondary electron trapping 

Trapping of electrons in electrostatic waves produces BGK modes which eventually 
collapse via the sideband instability. The sideband instability is due to the nonlinear 
oscillations of the electrons in the potential of the ESW. For electrons close to the 
bottom of the wave potential, their oscillation is that of a harmonic oscillator. The 
monochromatic bouncing frequency is Doppler shifted, due to the phase speed of the 
ESW. These sidebands can couple the electron energy to secondary high-frequency 
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Fig. 5: Contour plot of log(/e > 10~^) at time t = 200 for a system with a relativistic 
initial beam velocity of Vb = 0.9c. The exponential growth of the ESW is in the process 
of saturating via the trapping of electrons, forming a pair of BGK modes at p ~ ±3. 
The reduced significance of the higher temperature of species 4, with respect to species 
3, at high beam speeds allows the growth of forwards and backwards propagating ESWs 
simultaneously. Here x is given in units of 27r//c„. 

ESWs, which must have a wave number k < ku (Krasovsky |1994j ). The sideband 
instabiUty is a Umiting factor for the Ufetime of the ESW which has implications in 
the presence of an external magnetic field in particular. To demonstrate this we show 
in Fig. El the amplitude of the ESW driven by beam 1 at k = ku for beam speeds 
ranging from 0.15c to 0.9c. We find that the ESWs moving at nonrelativistic phase 
speeds saturate smoothly, which is in line with the wave saturation of the ESW in the 
Vlasov simulation in Dieckmann et al. ('2004b]). Their lifetime is significant and, 
in the presence of a weak external magnetic field B orthogonal to the wave vector 
k, the trapped electrons would undergo substantial ESA (Eliasson, Dieckmann Sz 
Shukla |2005j ). The ESA is proportional to Vph\B\ and the comparatively low Vph 
would limit the maximum energy the electrons can reach. 

The ESW driven by the beam with Vb = 0.6c grows to a larger amplitude and 
then collapses abruptly, which confirms the finding in Dieckmann et al. ( j2004bj l 
that the BGK modes become more unstable the larger the ratio between beam speed 
to electron thermal speed becomes. The Lorentz force eVpfi\B\ excerted by a BJ-k 
would be significantly higher than for the case Vh = 0.15, however the short lifetime 
of the wave would here prevent electrons from reaching highly relativistic speeds. 
By increasing the beam speed to Vb = 0.9c we obtain a stabilisation of the ESW, in 
line with the results in Dieckmann et al. ( j2()()4aj ). Here the Lorentz force is strong 
and the lifetime of the saturated ESW is long constituting a formidable electron 
accelerator, provided the field evolution is not strongly influenced by B. 

We now summarise the principal results from each beam speed. Note that in 
particular the ESWs driven by the mildly relativistic proton beams show oscillations 
after their initial saturation. To identify the origin of these fluctuations we apply a 
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Fig. 6: The time evolution of the ESW amphtudes with k = that are driven by beam 
1 for beam speeds of 0.15c, 0.2c, 0.6c and 0.9c. ESWs generated at non-relativistic beam 
speeds saturate smoothly and stabilise briefly. At Vb = 0.6 the ESW collapses abruptly 
and at relativistic beam speeds, Vb = 0.9c we observe a stabilisation of the ESW for a 
period of almost 1000ti;~g^, equivalent to almost 160 plasma periods. Approximate regions 
of stability are highlighted on each plot. 

Window Fourier Transform to the amplitudes of the ESWs with the wavenumbers 
ku/2 and 

3.3 Beam speed, vi, = 0.15c. 

In Fig. dwe find a strongly asymmetric ESW growth for Vb = 0.15c. Here the ther- 
mal spread of the individual plasma species is not small compared to Vb- Therefore 
the ESW driven by the cooler beam 1 grows and saturates first. It rapidly collapses 
and this collapse inhibits a further growth of the ESW driven by beam 2. The 
initial monochromatic ESW collapses into a broad wave continuum centred around 
CO = Wp.i. We find the equivalent broad wave continuum for k = ku/2 centred at 
Lj = These waves propagate at phase speeds comparable to Vb and represent 

the structures remaining of the initial BGK mode. The ESW spectrum at k = ku/2 
further shows two wave bands that are separated by a frequency of ujp^i from the 
main peak. These two wave bands could be pumped by a beat between the turbu- 
lent structure centred at a; = ujp^i/2, k = ku/2, the turbulent structure centred at 
Lu = Up^i and k = ku and the Langmuir wave with k = ku/2. Evidence for this is 
the correlation between both turbulent structures and the wave bands at t ~ 2500. 
At this time most wave power at a; ~ t^p,i/2 is absorbed at A; ~ ku/2. At the same 
time wave power at w ~ ujp^i is absorbed at k = ku while the power in the wave 
bands with uj « ujp^i/2 and to « 1.5a;p^i grows at k = ku/2. The faster of these two 
has a phase speed of co/k ^ 3vb = 0.45c. 

This beat wave is thus considerably faster than the initial ESW. By its large 
amplitude it could trap electrons. This is confirmed by Fig. |H1 where we find a 
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Fig. 7: The ESW spectrum for a beam speed of Vb = 0.15c: (a) shows the frequency 
power spectrum of kc = ku and (b) shows the frequency power spectrum of kc = ku/^ as 
a function of time. The colour scale shows log^g W^2('^5 1, kc)- The ESW driven by beam 1 
with Vb > corresponds to the wave with u = Up^i in (a). 

BGK mode in the electron distribution centred around a momentum p/nieC ~ 0.52 
which corresponds to a speed 0.46c. The fastest electrons of this BGK mode reach 
p/rrieC 0.68 or a speed of 0.56c 3.75^;,. 

3.4 Beam speed, vi, = 0.2c. 

The simulation with Vb = 0.2c shows a similar wave coupling as shown in Fig. 121 
We find a turbulent ESW spectrum that has been driven by the beam 1 and that 
covers waves with phase speeds centred at Vb with a spread that is a significant 
fraction of the beam speed. Initially an ESW is also driven by beam 2 at cj ~ ujp^i 
but, in line with Fig. [71 it collapses simultaneously with the ESW driven by beam 
1. The turbulent wave spectrum with oj ~ Up^i appears to couple with 

the equivalent spectrum at w ~ at k = ku/2 and the Langmuir wave with 

LJ = ±ujp^i and k = ku/2 to give wave bands at a; — 0.5u;p^i and at a; 1.5u;p^i. 
The phase speed of the ESW band at u; « 1.5u;p^i is to/k ^ 3vb and that of the band 
at w uJp^i/2 is uj/k Vb- 

The large amplitudes of both bands suggests that they might also be trapping 
electrons. This is confirmed by Fig. where we show the electron momentum 
distribution at t = 2000. We find BGK modes centred at p/rrieC « 0.75 or a speed 
of 0.6c. The BGK modes extend up to a peak momentum of p/mgC ~ 1 or a speed 
of 0.87c respectively equivalent to A.bvb- 

3.5 Beam speed, Vb = 0.4c. 

As we increase the beam speed to Vb = 0.4c the ESW driven by beam 1 stabilises 
and the collapsing ESW driven by beam 2 drives the broadband turbulence as 
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Fig. 8: Contour plot of log(/e > 10"^) for = 0.15c at the time t = 3500: The initial 
electron BGK modes have collapsed and they can be shown to form a plateau distribution 
at p/rrieC < 0.3. Centred at p/rrieC ~ 0.52 we find BGK modes driven by the beat wave. 
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Fig. 9: The ESW spectrum for a beam speed of Vb = 0.2c: (a) shows the frequency power 
spectrum of kc = fc„ and (b) shows the frequency power spectrum of kc = fc„/2 as a 
function of time. The colour scale shows \ogiQ W2{u!, t, kc). The ESW driven by beam 1 
with Vb > corresponds to the wave with u = Up^i in (a). 
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Fig. 10: Contour plot of log(/e > IQ-^) for Vb = 0.2c at the time t = 2000: The 
initial electron BGK modes have collapsed and form a plateau distribution. Centred at 
p/rrieC ~ 0.75 we find BGK modes driven by the beat wave. 



we see from Fig. 1111 The peak energy of the two counter-propagating ESWs is 
now comparable. Since we keep the thermal spread of the beams constant while 
increasing Vb, the thermal effects are reduced and the growth rate of both waves 
approaches the peak growth rate for the cold beam instability. We find the growth 
of ESWs at k = ku with ~ 2ujp^i. Since these ESWs have the same wave 
number as the initial wave they can not be produced by self-interaction of the 
initial ESW since this would also double the wave number. Instead we believe that 
these high frequency waves correspond to the Doppler shifted bouncing frequency of 
the trapped electrons in the potential of the strong ESW. The fastest electrons reach 
approximately twice the phase speed Vph of the ESW as can be seen, for example, 
at the lower beam speed of Vb = 0.2c in Fig. 0] These electrons are thus, due to 
the fixed wave number k^ of the BGK modes, interacting with secondary ESWs 
with CO ~ 2ujp^i. The ESWs with ~ 1.2ujp_i at k = ku/2 are, on the other hand, 
possibly due to a beat between the ESWs with |ti;| = ujp^i at k = ku and the ESWs 
with |a;| = at k = ku/2 and Langmuir waves, as for the slower ESWs. 



3.6 Beam speed, Vb = 0.6c. 

The ESW wave spectrum for Vb = 0.6c is similar to that for Vb = 0.4c. Again we 
find that both beams grow strong ESWs close to \uju\ = ^p,i but that it is the still 
stronger ESW driven by beam 1 that stabilises while that driven by beam 2 collapses 
into a broadband spectrum. This is also reflected by the ESW spectrum at fe„/2 
where we find the strongest wave activity at u; ~ —ujp^i/2. Two additional wave 
bands with frequencies u « tOp^i are observed at ku/2, i.e. long Langmuir waves 
are produced by the nonlinear processes. These waves have twice the phase speed 
of the initial ESWs and by their superluminal phase speed they can not interact 
resonantly with the electrons. At k = ku, on the other hand, we find high-frequency 
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Fig. 11: The ESW spectrum for a beam speed of Vb = 0.4c: (a) shows the frequency 
power spectrum of kc = and (b) shows the frequency power spectrum of kc = ku/2 as 
a function of time. The colour scale shows log^Q W2{uj, t, kc). The ESW driven by beam 1 
with Vb > corresponds to the wave with u = Up^i in (a). 



modes at \uj\ ^ l.Owp^i with a subluminal phase speed to/ku « 0.95c. 

The phase speeds of these ESWs is higher than the peak speed the electrons reach 
in the inital BGK mode as shown in Fig. 1131 We can thus not explain its growth 
by a streaming instability between the trapped electrons and, for example, the 
untrapped electrons. We can obtain phase speeds of the ESW bands that are higher 
than the speed of the trapped electron beam, however, by applying the relativistic 
Doppler shift to the electron bouncing frequency in the ESW wave potential. The 
rest frame of the ESWs moves with the speed Vph = 0.6c. The frequency of the 
ESWs in the observer frame is, according to Fig. I12L lOq ~ l.Qup^i. With the 
relativistic Doppler equation we would obtain a bouncing frequency of the electrons 

1 /2 

in the ESW frame of reference of iOb = ([1 — Vph/c]/[l + Vph/c]) ' uJq = O.Sup^i. We 
use the nonrelativistic estimate of the electron bouncing frequency in a parabolic 
electrostatic potential = ekuE/nie and the corresponding width of the trapped 
electron island vf,^ = 2eE /nieku to eliminate the electric field E. We obtain the 
relation oJb/ku = vtr/V^- Since for our cold plasma species the velocity width of the 
island of trapped electrons must be comparable to Vph to trap the bulk electrons we 
get an estimate for ojb ~ u!p^i/\/2 which is close to ujo- We may thus indeed interpret 
the two sidebands observed in Fig. ^Jas the Doppler shifted bouncing frequency of 
the electrons. 

These sideband modes driven by the beams of trapped electrons have a phase 
speed that is just below c and a large amplitude. Since their phase speed is compara- 
ble to the fastest speed the electrons reach in the initial BGK mode, they should be 
capable of trapping some of these electrons. This is confirmed by Fig. 1141 where we 
find BGK modes centred at the momentum p/nieC ~ 3.4 which accelerate electrons 
up to the peak momentum p/nieC 7 or a speed of 0.99c. 
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Fig. 12: The ESW spectrum for a beam speed of Vb = 0.6c: (a) shows the frequency 
power spectrum of kc = ku and (b) shows the frequency power spectrum of kc = ku/2 as 
a function of time. The colour scale shows log^g W2{uj, t, kc). The ESW driven by beam 1 
with Vb > corresponds to the wave with u = Up^i in (a). 




Fig. 13: Contour plot of log(/e > 10^^) for Vb = 0.6c at the (early) time t = 200: The 
initial electron BGK modes have just been formed and the trapped electrons reach a 
momentum up to p/rrigC ~ 3. 
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Fig. 14: The electron momentum distribution for vi, = 0.6c at the (late) time t = 3500: 
The initial electron BGK modes have collapsed and the electrons that are trapped by the 
sideband mode reach a momentum up to p/rrieC ~ 7. 

3.7 Beam speed, Vb = 0.8c and Vb = 0.9c. 

As we increase the beam speed to = 0.8c the qualitative evolution of the ESWs 
changes. At this high beam speed the thermal spread of the plasma species is 
negligible. The growth rates of the waves for both beams is approximately that of 
the cold beam instability and the ESWs grow symmetrically. Both saturated ESWs 
are stable. After the saturation each ESW shows a sideband, similar to that in Fig. 
1131 at a frequency modulus lo ~ l.Swp. The phase speeds of these sideband modes 
are presumably just above c. We find a growing second mode at A; = ku/'2 at the 
frequency modulus uj k, O.Swp^i which moves at a superluminal phase speed. 

We observe the same growth of sideband modes for the fastest beam speed of 
V}j = 0.9c in Fig. ^1 The sideband modes at k = ku have a frequency modulus 
Lj PS 1.2a;p^i and phase speeds just above c. As for Vb = 0.8c the sideband mode at 
k = ku/2 has a frequency of a; ~ O.Swp^i and thus a superluminal phase speed. 

For both beam speeds Vb = 0.8c and Vb = 0.9c the sideband modes appear to 
have a superluminal phase speed and they can therefore not trap electrons. No 
secondary BGK modes should develop for these beam speeds. This is confirmed by 
Fig. 1171 where we show the phase space distribution of the electrons for Vb = 0.9c 
at the simulation's end. The phase space distribution at high momenta shows no 
evidence of a BGK mode despite the strong sideband modes in Fig. 1161 

In contrast to the PIC simulations in Dieckmann et al. ( (2004aj ) the Vlasov 
simulation code shows the growth of sideband modes and what appears to be waves 
resulting from a parametric instability. These secondary waves grow to a large am- 
plitude at which they can nonlinearly interact with the electrons. For the beam 
speeds up to Vb = 0.4c the waves generated by the parametric interaction have 
been strongest. Here the turbulent wave fields interact with the waves with ujp^i 
to produce a wave with a higher frequency. In contrast to plasma beat wave ac- 
celerators, which have recently been reviewed by Bingham, Mendonca & Shukla 
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Fig. 15: The ESW spectrum for a beam speed of Vb = 0.8c: (a) shows the frequency 
power spectrum of kc = and (b) shows the frequency power spectrum of kc = ku/^ as 
a function of time. The colour scale shows log^Q W2{uj, t, kc)- The ESW driven by beam 1 
with Vb > corresponds to the wave with uj = ujp^i in (a). 
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Fig. 16: The ESW spectrum for a beam speed of Vb = 0.9c: (a) shows the frequency 
power spectrum of kc = ku and (b) shows the frequency power spectrum of kc = ku/2 as 
a function of time. The colour scale shows log^g W^2('^, t, kc). The ESW driven by beam 1 
with Vb > corresponds to the wave with u = Up^i in (b). 
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Fig. 17: Contour plot of log(/e > 10"^) for Vb = 0.9c at the time t = 4000: No well-defined 
BGK mode is visible in the electron distribution since the sideband unstable modes have 
a superluminal phase speed. 

f |2UL)4j '). for which two high-frequency electromagnetic waves beat to yield a low 
frequency ESW that can accelerate the electrons, our parametric coupling couples 
low frequency ESWs to an ESW with a higher frequency. Its larger phase speed 
can accelerate the trapped electrons to higher peak speeds. We may thus call it the 
"inverse plasma beat wave accelerator". For a beam speed of = 0.6c the turbulent 
wave fields do not noticably interact with ujp^i. Instead a sideband unstable mode 
with k = ku develops, i.e. at the largest allowed wave number for nonrelativistic 
BGK modes (Krasovsky |1994j ). This mode has for Vh < 0.6c a phase speed below 
c and the electron phase space distribution shows a fast BGK mode driven by it. 
For even higher Vj, the probably superluminal phase speed of the sideband modes 
suppresses their interaction with the electrons. 

The complex spectrum of secondary waves and their nonlinear interactions with 
the electrons, which has not been observed clearly by Dieckmann et al. ( i20 0laj ) . 
suggests a stronger dependence of the electron heating on than for the simulations 
by Dieckmann et al. ( |2004aj ) . There the momentum distributions could be matched 
if one were to scale the momentum axis to Vhj{vi,). We thus do the same here and 
we integrate the electron phase space distribution over x. The result is shown in 
Fig. d 

We find that the final momentum distributions in the chosen normalization of 
the p-axis agree well up to Vf, = 0.6c and for densities larger than 10~^. Here the 
Vlasov code results are similar to the equivalent PIC simulations in Dieckmann et 
al. ( j2()()4aj ). We find, however, significant differences at lower densities, which are 
not well-represented by the PIC code, and at the corresponding higher momenta. 
The parametric instabilites and the sideband instabilities have further accelerated 
electrons beyond the peak momenta measured in Dieckmann et al. f j2004aj '). This 
is particularly pronounced for the beam speed Vb = 0.6 in which we find an electron 
density plateau extending to a value of -|-10 for the normalised momentum, i. e. 



102 



10° 



10-2 



10-4 



10-6 

-15 -10 -5 5 10 15 

Fig. 18: The electron momentum distributions at the simulation's end times: 

twice as high as the corresponding value at negative momenta. For these beam 
speeds the secondary waves have been most efficient as electron accelerators. A 
further increase of beyond 0.6c yields broadening momentum distributions. The 
peak momentum in units of rUeC is, for positive momenta comparable for Vb = 0.6c 
and for Vb = 0.8c. The peak momentum for Vb = 0.9c is about twice as high as for 
Vb = 0.6c. 

The peak relativistic kinetic energies K = mc^e~^(7 — 1) in eV the electrons 
reach are K{vb = 0.15c) = 1.3 x 10^ eV, K{vb = 0.2c) = 2.2 x 10^ eV, K{vb = 
0.4c) = 5 X 10^ eV, K{vb = 0.6c) = 3.2 x 10^ eV, K{vb = 0.8c) = 3.6 x 10^ eV 
and K{vb = 0.9c) = 8 x 10^ eV. Note that all these peak electron energies are 
comparable or above the threshold energy of 10^ eV: the injection energy for Fermi 
acceleration at perpendicular shocks, given by Treumann & Terasawa f |2001| ). 

4 Discussion 

The observation of the emission of highly energetic cosmic ray particles by SNRs 
suggests the acceleration of particles from the thermal pool of the ISM plasma 
to highly relativistic energies by such objects. The acceleration site is apparently 
linked to the shock that develops as the supernova blast shell encounters the ambient 
plasma (Lazendic |2004) ). Such shocks are believed to accelerate electrons and ions 
to highly relativistic energies by means of Fermi acceleration (Fermi jl949j : Fermi 
|19B4j ). The Fermi acceleration of electrons is most efficient if the shock is quasi- 
perpendicular (Galeev |1984j ). For such shocks, however, Fermi acceleration works 
only if we find a relativistically hot electron population prior to the shock encounter, 
since slow electrons could not repeatedly cross the shock and pick up energy. Since 
the plasma, into which the SNR shock expands, has a thermal speed comparable 
to that of the ISM or the stellar wind of the progenitor star with temperatures 
of up to a few eV, if we take the solar wind as reference, no mildly relativistic 
electrons may exist. A mechanism is thus required that accelerates electrons up to 




speeds at which their Larmor radius exceeds the shock thickness. As Galeev ( jl984| ) 
proposed, electrons could be pre-accelerated from the initial thermal pool to mildly 
relativistic energies by their interaction with strong ESWs in the foreshock region 
which, in turn, are driven by shock-reflected beams of ions. 

We have examined in this work the growth, saturation and collapse of ESWs 
in a system dominated by the presence of two counter-propagating proton beams. 
The currents of both beams cancel, allowing the introduction of periodic boundary 
conditions. These initial conditions can be motivated as follows: The upstream 
protons, that have initially been reflected by the shock, are rotated by the global 
magnetic field oriented perpendicularly to the shock normal and return as a second 
counter propagating proton beam, with a slightly increased temperature. The sys- 
tem modelled in this work represents a small region ahead of a SNR shock. Here 
the magnetic field has been neglected in order to focus on the ESWs and nonlin- 
ear BGK modes. While this does not describe a complete model for the foreshock 
dynamics of high Mach number shocks, it is applicable to parts of the foreshock 
of perpendicular shocks where local variations magnetic field cause it to vanish, or 
become beam-aligned. These variations may, for example, be due to the presence 
of turbulent magnetic field structures in the foreshock of SNR shocks (Jun & Jones 
[TM?] : Lazendic et al. 2004 J . 

The system is unstable to the relativistic Buneman instability (Thode & Sudan 
jl97,Sj ) which saturates via the trapping of electrons to form BGK modes (Rosen- 
zweig jl988| '). These trapped particle distributions are themselves unstable to the 
sideband instability and collapse after a period of stability. 

Since our model assumes that the particle trajectories are not affected by the 
local magnetic field and because we do not consider here electromagnetic waves or 
waves in magnetised plasma, we are able to utilise a relativistic, electrostatic Vlasov 
code. Previous work, for example (Thode & Sudan |1973j : Dieckmann et al. |2000j : 
Dieckmann et al. |2004bj : Shimada &: Hoshino |[2DD^ ) has made extensive use of PIC 
codes for this problem and in particular Dieckmann et al. ( |20n4bj ) have compared 
Vlasov and PIC codes in certain conditions. We benefit from the Eulerian Vlasov 
code's ability to resolve electron and ion phase space accurately, irrespective of the 
local particle density. This allows us to identify a secondary acceleration mechanism 
which may not be immediately apparent otherwise. Overall the results of this work 
are in agreement with previous studies (Dieckmann et al. J2000'; Dieckmann et al. 
|2004bj ) showing the lifetimes of the BGK modes to be dependent on the initial beam 
velocity. As Vb is increased, we observe a reduced ESW stability up to Vb = 0.8c but 
with a significantly increased stability at the highest beam speed. 

At low beam velocities {vb = 0.15c, 0.2c, 0.4c, 0.6c) we observe long wavelength, 
high phase- velocity modes. These are able to trap electrons, producing a population 
with kinetic energies above the injection energy for Fermi acceleration, even at non- 
relativistic beam, and thus shock, velocities. We believe these secondary (that is to 
say, not associated with the initial ESW saturation) trapped distributions to be the 
result of trapping in ESWs produced by parametric coupling between low frequency 
oscillations and plasma waves at Above Vb = 0.6c the ESWs produced by 

this coupling have super-luminal phase velocities and are unable to trap electrons. 
Hence we do not observe such BGK modes in the case of Vb = 0.8c, or 0.9c. Our 
simulation box, at 4:Tr /ku in length, can only accommodate one wavemode with 
wavenumber below that of the most unstable mode and this may have an infiuence 
on the appearance of this secondary acceleration. 

Future work has to examine how these parametric instabilities depend on a 



magnetic field and on the introduction of a second spatial dimension. We will need 
to consider significantly larger simulation boxes, capable of resolving a broader 
spectrum below ku- It may be the case that the availability of modes with k < ku 
will result in the partition of ESW energy across a greater region of the spectrum, 
perhaps inhibiting the trapping of electrons at high velocity. However, it may be the 
case that our observed coupling and resultant electron trapping is the first step of a 
cascade, capable of accelerating electrons to high energies for relatively modest shock 
velocities. This is since, even for Vh as low as 0.15c or corresponding shock speeds of 
7.5 X 10~^c which can be reached by the fastest SNR main shocks (Kulkarni et al. 
jl998p . electrons can reach energies of 10^ eV. According to Treumann & Terasawa 
( |2()()lj ). this may increase the electron gyroradius beyond the shock thickness by 
which they can repeatedly cross the shock front. These repeated shock crossings 
allow the electrons to undergo Fermi acceleration to highly relativistic speeds. Even 
higher energies could be achieved if we were to get shock precursors that outran the 
main shock as has been observed for the supernova SN1998bw (Kulkarni et al. 
jl998j ^. The numerical simulations in this work thus present strong evidence for 
the ability of ESWs and processes driven by electrostatic turbulence to accelerate 
electrons beyond the threshold energy at which they can undergo Fermi acceleration 
as it has previously been proposed, for example by Galeev ( |1984j ). 
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